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Assuming that the final state ol hadronization takes place along the freezeout line, which is defined 
by a constant entropy density, the antiproton-to-proton ratios produced in heavy-ion collisions are 
studied in framework of the hadron resonance gas (HRG) model. A phase transition from quark- 
gluon plasma to hadrons, a hadronization, has been conjectured in order to allow modifications 
in the phase space volume and thus in single-particle distribution function. Implementing both 
modifications in the grand-canonical partition function and taking into account the experimental 
acceptance in heavy-ion collisions, the antiproton-to-proton ratios over center-of-mass energies y/s 
^ ' ranging from AGS to RHIC are very well reproduced by the HRG model. Comparing with the 

same particle ratios in pp collisions results in a gradually narrowing discrepancy with increasing 
y/s. At LHC energy, the ALICE antiproton-to-proton ratios in pp collisions turn to be very well 
described by HRG model as well. It is likely that the ALICE heavy-ion program will produce 
\ the same antiproton-to-proton ratios as the pp program. Furthermore, the ratio gets very close to 

C^l . unity indicating that the matter-antimatter asymmetry nearly vanishes. The chemical potential 

calculated at this energy strengthens the assumption of almost fully matter-antimatter symmetry 
at LHC energy. 
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I. INTRODUCTION 



> . 

£N| , The collective properties of hot and dense QCD matter is one of the main objectives of the heavy-ion program. 
The possible modification of the transport coefficients, like phase structure and effective degrees of freedom, with 
increasing incident energy and system size likely provides fruitful tools to study the collective properties. Furthermore, 
it comprehensively characterizes the particle production itself [l[ . The universal description for the particle production 
at different center-of-mass energies of nucleon-nucleon collision y/s can be achieved by studying the fluctuations [l|, 0] 
and multiplicities of particle ratios [H, H|, which are not only able to determine the freezeout parameters, temperature 
T and chemical potential //, but also - to a large extend - eliminate the volume fluctuations and the dependence of the 
freezeout surface on the initial conditions. The fluctuations in the particle production are also essential regarding with 
examining the existing statistical models 0, Q , characterizing the particle equilibration Q and search for unambiguous 
k>< i signals for the creation of the quark-gluon plasma (QGP) 

Combining various fluctuations and multiplicities of integrated particle ratios [H, H[ makes it possible to determine 
d ' the freezeout parameters by means of thermal models. The search for common properties of the freezeout parameters 
in heavy-ion collisions at different center-of-mass energies has a long tradition [10|, Hll. Different models have been 
suggested [lOrll^ . In present work, the constant entropy density s normalized to T 3 jl2| is used. It implies that the 
hadronic matter in rest frame of produced system is characterized by constant degrees of freedom. The scaling of the 
quantity s/T 3 by 7r 2 /4, which reflects the phase space volume, obviously results in the effective degrees of freedom of 
free gas [Fa ]. 

On the other hand, it has been concluded that none of statistical models [l|| would be able to perfectly describe the 
experimental data. Using grand-canonical ensembles and phase space modification factor Q leads to a much better 
description of the experimental data [13, EH ■ The hadron resonance gas (HRG) model turns to fairly reproduce the 
particle multiplicities, the dynamical fluctuations of various particle ratios [ijj and the thermodynamics of the strongly 
interacting system. In doing this, Q plays an essential role. It has been noticed that the fluctuations over the whole 
range of y/s exhibit a non-monotonic behavior and Q is varying with increasing y/s (l7l . [Lsj . 

Enhanced antibaryons are conjectured as indicators for the formation of deconnned QGP [2(| , whereas the possible 
annihilation might suppress such enhancement pll ] . The collision initial conditions and formation time can be reflected 
by surviving antibaryons, where antiproton is a good candidate. Decelerating or even stopping of incident particle and 
its break up in inelastic collisions have been discussed in literature (22j . Therefore, np/n p ratios can be used to study 
the baryon transport and production and therefore would have significant astrophysical consequences. By mentioning 
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astrophysics, it seems in order to recall now that the np/n p ratios have been calculated and also observed in the cosmic 
rays revealing essential details on astrophysical phenomena [23[ . In present work, n p / n p ratios are calculated in HRG 
model and compared with the heavy-ion collisions at y/s ranging from AGS to LHC. Thus, np/n p ratios and their 
description in the thermal model seem to provide fruitful tools to study the evolution of matter-antimatter asymmetry 
with changing energy. In section HH the model is introduced. Section UTTl is devoted to the results and conclusions. 



II. HADRON RESONANCE GAS MODEL 



The hadron resonances treated as a free gas [6, 24 27] arc conjectured to add to the thermodynamic quantities in the 
hadronic phase of heavy-ion collisions. This statement is valid for free as well as strong interactions between the hadron 
resonances themselves. It has been shown that the thermodynamics of strongly interacting system can be approximated 
to an ideal gas composed of hadron resonances [26l - [28| . Also, it has been shown that the thermodynamics of strongly 
interacting system can be approximated to an ideal gas composed of hadron resonances with masses < 1.8 GeV [26l — 
|28| . The heavier constituents, the smaller thermodynamic quantities. The main motivation of using the Hamiltonian 
is that it contains all relevant degrees of freedom of confined and strongly interacting matter. It implicitly includes 
the interactions that result in the formation of new resonances. In addition, HRG model is used to provide a quite 
satisfactory description of particle production and collective properties in heavy-ion collisions [13, fl7Hl9l . [24T - (28j . 

At finite temperature T, strange /is and isospin fij a and baryochemical potential /is, the partition function of one 
single hadron reads [Tt], fig 



In Z gc (T, V, /i) = ±^2 V J edk In I 1 ± 
lnZ gc (T,V,/i) = 5Z±^ V J k 2 dk In M ± 7 Q exp 



1+ l^p^-a 



l/(9-l) N 



/i. t - ej(k) 
T 



(1) 
(2) 



where £j(fc) = (fc 2 + m 2 ) 1 ^ 2 is the i-th particle's dispersion relation (single-particle energy) in equilibrium and ± stands 
for bosons and fermions, respectively, g is the spin-isospin degeneracy factor and 7 = 7^7™ stand for the quark phase 
space occupancy parameters, where n and m being number of light and strange quarks, respectively. A = exp(c/x/T) 
is the fugacity, where \x is the chemical potential multiplied by corresponding charge c, like baryon, light and strange 
quark numbers, etc. Summing over all hadron resonances results in the final thermodynamics in the hadronic phase, 
since no phase transition is conjectured in HRG model. 

Parameters a and /3, are Lagrange multipliers in the entropy maximization. The physical meaning of a is a controller 
over the number of particles in the phase space, i.e, acting as chemical potential, a = In £ — In T — In N . It also combines 
intensive variables, T and N with an extensive one £ = gi^i 1 where q is energy density of i-th cell in the phase 
space of interest. The most probable state density can be found by the Lagrange multipliers, where one of them, a, 
has been expressed in term of the second one, /3, and the occupation numbers of the system. Apparently, a gives how 
the energy £ is distributed in the microstates of the equilibrium system and therefore, can be understood as another 
factor controlling the number of occupied states at microcanonical level. 

The generic parameter Q(x, k) is conjectured to reflect the change in phase-space, when the hadronic degrees of 
freedom are replaced by partonic ones and vice versa. Also, it can be interpreted as a measure for the non-extensivity. 
Comments on the statistical parameter 7 = exp(— a) are now in order. It gives the averaged occupancy of the phase 
space relative to equilibrium limit. Assuming finite time evolution of the system, 7.; can be seen as a ratio of the 
change in particle number before and after the chemical freeze-out, i.e. 7, = m(t)/rii(oo). The chemical frcczeout is 
defined as a time scale, at which there is no longer particle production and the collisions become entirely elastic. In 
case of phase transition, 7 j is expected to be larger than one, because of the large degrees of freedom, weak coupling 
and expanding phase space above the critical temperature to QGP. 

The quark chemistry is given by relating the hadronic chemical potentials and 7 to the quark constituents. 7 = 7g 7™ 
with n and m being the number of light and strange quarks, respectively. [Ib = 3/i g and [is = n q — n s , with q and 
s are the light and strange quark quantum number, respectively. The baryochemical potential for the light quarks is 
\i q = (p u + /i(j)/2. The iso-spin chemical potential /(/ 3 = (p, u — /id)/2. Therefore, vanishing /i/ 3 obviously means that 
light quarks are conjectured to be degenerate and fi q = fi u = /id- The strangeness chemical potential, fig, is calculated 
in HRG model as a function of T and (13- At finite T and hb, (J-s is calculated under the condition of the overall 
strangeness conservation. 

Absolving chemical freezeout means that the hadrons finally decay to stable hadrons or other resonances 
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where bj^ being decay branching ratio of j'-th hadron resonance into z-th stable particle. The particle number density 
in grand-canonical ensemble is no longer constant. In an ensemble of N r hadron resonances. 



E 



2tt 2 



k 2 dk- 



exp 



(4) 



In carrying out the calculations given in Fig. [TJ a full grand-canonical statistical set of the thermodynamic parameters 
has been used. Corrections due to van der Waals repulsive interactions has not been taken into account [26| . Seeking 
for simplicity, one can for a moment assume the Boltzmann approximations, especially at ultra-relativistic energy, in 
order find out the parameters, on which np/n p ratio depends. Then, at finite isospin fugacity A/ 3 , 



(^1 X d ) 



(5) 



III. RESULTS AND CONCLUSIONS 



As mentioned above, the formation of deconfincd QGP strongly depends on the initial conditions in the early stage 
of heavy-ion collisions, where the baryon number transport (deceleration of the incoming particle) or even its entire 
stopping is assumed to be achieved [29j. Therefore, it is widely conjectured that such deconfined matter would not 
be produced, especially in pp collisions at low energies. On the other hand, the degrees of baryon number transport 
likely affects the overall dynamical evolution of the system. The initial parton equilibrium [30(, thermal/chemical 
equilibrium (3l| . time evolution [32j and the final particle production [33| are affected as well. 

In Fig. [T] n p /n p ratios are given in dependence on center-of-mass-energy -y/s ranging from AGS to LHC. The pp 
collisions (empty symbols) are compared with the heavy- ion collisions (solid symbols). In the nucleus-nucleus collisions, 
the spatial and time evolution of the system is too short to assure initial conditions required in order to derive the 
parton matter into deconfined QGP. We notice that np/n p ratios up to RHIC energy are smaller than unity indicating 
an overall excess of p over p. This would also imply that a certain fraction of the baryon number is likely transported 
from the incident projectile at beam rapidity to mid-rapidity region of the system. Thus up to RHIC energy, the 
mid-rapidity region seems to have finite net-baryon number. 

As mentioned above, HRG model has been successfully utilized to describe heavy-ion collisions and the thermody- 
namics of lattice QCD as well 0, [T3, EH [H|-[26j]. The collective flow of the matter in heavy- ion collisions makes it 
obvious that HRG model underestimates the np/n p ratios measured in the pp collisions, especially at low energies. 
On the other hand, the ratios from the heavy-ion collisions [34[ are very well reproduced by HRG model. In both 
types of collisions, there is a dramatic increase in the np/n p ratio. For instance, in RHIC heavy- ion collisions, the 
ratio is 0.65 ± 0.05 U/s = 130 GeV), which is obviously much higher than the ones measured at SPS (~ 0.07 ± 10% 
at y/s = 17 GeV) [aS] and AGS (~ 2.5 x 10" 4 ± 10% at y/s = 5 GeV) [H]. Also in the pp collisions, the n p /n p ratios 
obviously increase from NA49 over ISR and 200 GcV-RHIC to ALICE. Up to RHIC energy, it is obvious that the 
matter-antimatter asymmetry decreases from nearly 100% at AGS to 65% and 130GeV-RHIC. At LHC energy, while 
the differences between n p /n p ratios of pp and heavy- ion collisions apparently disappear, their values approach unity, 
i.e. almost 0% matter-antimatter asymmetry. 

The dashed line in Fig. [T]gives the fitting according to Regge model, which reasonably describe the n p /n p data from 



different pp collisions. This fitting is elaborated in Ref. [37] . The baryon pair production is assumed to be governed 



by pomeron exchange and baryon transport by string-junction exchange [38|. Therefore, it has been concluded that 
ALICE results seem to be consistent with the baryon number transport. 

The solid line gives the results from HRG model, section [TTJ No fitting has been utilized. Over the whole range 
of y/s, the two parameters 7 and Q remain constant, 1.0 and 0.5, respectively. It seems to be a kind of a universal 
curve. It describes very well the heavy-ion results. Also ALICE pp results are reproduced by means of HRG model. 
So far, it can be concluded that the future heavy- ion program of ALICE likely will produce n p /n p ratios similar, 
when not entirely identical, to the ones measured in this current pp program. Secondly, the ratios run very close to 
unity implying almost vanishing matter-antimatter asymmetry. Some astrophysical consequences of these results will 
shortly be elaborated below. On the other hand, it can also be concluded that the thermal models including HRG 
seem to be able to perfectly describe the hadronization at very large energies and the condition deriving the chemical 
freezeout at the final state of hadronization, the constant degrees of freedom or S(y/s, T) — 7(A/tt 2 )VT 3 , seems to be 
valid at all center-of-mass-energies raging from AGS to LHC. 

Expression (|5|) can be used to calculate the antiproton (antiquark) chemical potential. At a vanishing isospin 
chemical potential /i/ 3 , i. e. degenerate light quarks 

~ exp(— 2/Xp/Tcfc) = exp(-6 / ti ? /T c / l ), (6) 

71-n 
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Fig. 1: np/n p ratios depicted in the whole available range of ^/s. Open symbols stand for the results from various 
pp experiments (labeled). The solid symbols give the heavy-ion results from AGS, SPS and RHIC, respectively. The 



fitting of pp results according to Regge model is given by the dashed curve [37[ • The solid curve is the HRG results 
at Q = 0.5 and 7 = 1.0, which apparently perfectly describes all heavy-ion collisions, besides the ALICE pp results. 
Extending it to the LHC energy obviously shows that the ratio itself is very much close to unity. Contrary to the 
dashed curve, the solid line is not a fitting to experimental data. 



where T c h is the freezeout temperature. Substituting np/n p of ALICE experiment at 7 GeV and T c h — 0.175±0.02 GeV 
in Eq. © results is 

Hp = 3fiq ~ 0.537 MeV (7) 

The ratio of baryon density asymmetry to photon density, 77, has been measured in WMAP data [39| . then rip — n p 
from ALICE experiment at 7 GeV gives a photon number density ranging between 

4.4 x 10 4 < n 7 < 7.7 x 10 4 . (8) 

When modelling the volume of the early universe at energies equal to top ALICE, 7TeV, the photon number, _/V 7 = 
n-y V, can be calculated. It is obvious that the resulting iV 7 will be extraordinarily large indicating dominant radiation 
against all other equations of state in this era where matter and antimatter are supposed to be almost symmetric. 

The matter-antimatter asymmetry is one of the greatest mysteries in modern physics. The well-known Sakharov's 
conditions pol ]; baryon number violation and C and CP violation besides out-of-equilibrium processes; are conjectured 
to give a solid framework to explain why should matter become dominant against antimatter with the universe cooling- 
down of expansion. Recently, a fourth condition has been suggested; b-l violation, where b and I being baryon and lepton 
quantum number, respectively. In standard model, both b and I are not conserved, explicitly. Non-perturbatively, 
they can be violated through sphalerons, for instance, which derive the two-direction conversion b f-> I with a certain 
number of selection rules |41[ . That the heavy-ion program seems to give systematic tools to survey the evolution of 
matter-antimatter asymmetry from ~ 0% to ~ 100% in the nuclear interactions would make it possible to devote such 
experimental facilities to study the evolution of Sakharov's conditions. 
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